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We study brane-world preheating in massive chaotic inflationary scenario where scalar fields are 
confined on the 3-brane. Assuming that quadratic contribution in energy densities dominates the 
Hubble expansion rate during preheating, the amplitude of inflaton decreases slowly relative to the 
standard dust-dominated case. This leads to an efficient production of \ particles via nonpertur- 
bative decay of inflaton even if its coupling is of order g = 10 . We also discuss massive particle 
creation heavier than inflaton, which may play important roles for the baryo- and lepto-genesis 
scenarios. 
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I. INTRODUCTION 

Recent studies of superstring or M-theory predict that 
our universe may consist of a brane world embedded in a 
higher-dimensional spacetime In particular, it was 

shown that the 10-dimensional Eg x E$ heterotic string 
theory, which is a strong candidate to describe our real 
world, is equivalent to an 11-dimensional M-theory com- 
pactified to M 10 x S 1 /Z 2 f|. Then the 10-dimensional 
spacetime is expected to be compactified into M 4 x CY 6 . 
Therein the standard model particles are confined on the 
3-dimensional brane. On the other hand, gravitons prop- 
agate in the full spacetime. 

In conventional Kaluza-Klein theories, extra dimen- 
sions are compactified on some manifolds in order to ob- 
tain 4-dimensional effective gravity theories. In contrast, 
Randall and Sundrum Q recently proposed a new type 
of compactification, in which extra dimension is not com- 
pact but gravity is effectively 3-dimensional. Their model 
indicates that we are living on the 3-brane with positive 
tension embedded in 5-dimensional anti de-Sitter bulk 
spacetime. This idea invokes a particular interest for the 
evolution of the early universe, which may exhibit some 
deviations from the standard cosmology by the modifi- 
cations of 4-dimensional Einstein equations on the brane 

In this respect, cosmological inflation was investigated 
by several authors in the context of the brane-world 
scenario where scalar fields are confined on the 3-branc. 
The effect of the brane in higher dimensions induces ex- 
tra terms in 4-dimensional effective gravity especially for 
high energies, e.g., a quadratic term of energy density, 
which leads to the larger amount of inflation relative to 
standard inflationary scenarios. 

If the modifications due to extra dimensions persist in 
the post inflationary phase, this may nontrivially affect 
on the process of reheating. In the early stage of re- 
heating, particles coupled to inflaton are efficiently cre- 



ated by the nonperturbative and nonequilibrium process 
which is called preheating |9|-|l2[. Since the particle pro- 
duction during preheating is sensitive to the background 
evolution pCfl , it is expected that the standard picture 
of preheating will be altered due to the modifications in 
Einstein equations on the brane. In fact we will show 
that the brane-world preheating broadens the parame- 
ter ranges of coupling constants where particles are effi- 
ciently produced. 



II. BRANE INFLATION 

In the brane world scenario by Randall- Sundrum [Q, 
where gravity as well as matter field are effectively con- 
fined on the 3-brane in the 5-dimensional spacetime, the 
5-dimcnsional Einstein equations can be recast in the fol- 
lowing gravitational equations on the brane Q : 
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(2.1) 



where T^ v and Tt^ v represent the energy-momentum ten- 
sor on the brane and the quadratic term in T^ v , re- 
spectively. E^ is a part of 5-dimensional Weyl tensor, 
which carries the information in the bulk. Note that 
the four- and five-dimensional gravitational constants, 
Mpi and M5, are related with the 3-brane tension, A, as 
A = 487rMf /Mpj.f] Here the 4-dimensional cosmological 
constant A4 is assumed to be zero. 

Let us analyze the gravitational equation ( J2 . L[ ) in 
cosmological contexts. Adopting a flat Friedmann- 



*Our notation of the 5-dimensional Planck mass, M5, does 
not include the \/87r factor, which is different from the defi- 
nition in Ref. 
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Robertson- Walker (FRW) metric as a background space- 
time on the brane, we find the following Friedmann equa- 
tion §§ 



H 2 = 



8tt 
3M^' 



(2.2) 



where a and p are the scale factor and the energy density 
of the matter on the brane, respectively. We have ig- 
nored the so-called "dark" radiation E^, which rapidly 
decreases as ~ a~ 4 during inflation even if it exists. 

At high energies the p 2 term is expected to play im- 
portant roles to determine the evolution of the Universe. 
Hereafter we will assume that the quadratic term domi- 
nates during inflation and preheating, i.e., the condition, 
p > 2A, holds until the end of preheating. The infla- 
ton field confined on the brane satisfies the Klein-Gordon 
equation, 



■ZH(l> + V'(<t>) = 0. 



(2.3) 



Note tha t during inflation the quadratic contribution in 
Eq. (2.2) increases the Hubble expansion rate, which 
makes the evolution of inflaton slower by Eq. (2.3). An- 



other Einstein equation in the quadratic-term dominant 
stage yields 
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This implies that inflation ends around 2V ~ 5(f) 2 . In 
the massive chaotic inflationary scenario with a potential, 
V(4>) — ^m 2 4> 2 , this condition gives 
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where we used the slow-roll condition, 



-4M 5 3 /V> in Eq. 

The inflaton mass, to, should be constrained by com- 
paring the amplitude of density perturbations produced 
during inflation with the CMB anisotropy. This is typ- 
ically done by evaluating the curvature perturbation, 
around the 55 e-foldings before the end of inflation, if £ is 
conserved on large scales. It was shown in Ref. jl3| that 
the conservation of £ holds for adiabatic perturbations 
irrespective of the form of gravitational equations by con- 
sidering the local conservation of the energy-momentum 
tensor. In the presence of a second scalar field, Xi curva- 
ture perturbations include isocurvature (entropy) pertur- 
bations jl3 14|, which provides the source of nonadiabatic 
pressure perturbations. However, as long as x does not 
exhibit strong growth during inflation as in the nonmini- 
mally coupled case 15 1, £ remains conserved even in the 
multi- field case jl6|. In fact, in the massive chaotic infla- 
tionary scenario with the interaction \g 2 <\y L x' \ X 1S sup- 
pressed during inflation for large couplings g required for 
the x particle production during preheating [17]], which 



ensures the conservation of £ on large scales. Although 
we will use this result here, the Weyl tensor E^ v will bring 
the 5-dimensional information into the brane, which may 
alter the result flisfl . More careful analysis will be re- 
quired regarding the evolution of £ during inflation [fill . 
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FIG. 1. The evolution of the inflaton condensate, 
4> — (m/Aft) 1 ' 2 (ft/Ms, during inflation and reheating. We 
start integrating with the initial value i n E q. (2.6). Inflation 
ends around (j> — 3 as estimated by Eq. (2.5). The system en- 
ters the reheating stage around cj> ~ 1. Inset: The evolution 
of the number of e-foldings, N = J F Hdt. We find that 55 
e-foldings are achieved during inflation. 

Calculating the number of e-foldings: N = J t ' F Hdt, 
and setting TV — 55, we obtain the value of inflaton when 
scales relevant for the CMB anisotropy crossed outside 
the Hubble radius, 
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The spectra of density perturbations at horizon re- 
entry can be evaluated as 8h = {2 / b)\C,\k= a H = 
(2/5)H5(f>/4i\k=aH , where the subscript, k = aH, denotes 
the values at first horizon crossing. Making use of the re- 
lation S(f> ~ H/2ir at first horizon crossings, we find 
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The COBE normalization gives 5h 
constrains the mass of inflaton as 

to/Mr ~ 1.4 x 1CT 4 . 



(2.7) 



2 x 10 '°. which 



(2. 



The ratio of Eq. (2.8) is crucial to determine the effi- 



ciency of parametric resonance during preheating. We 



' While the formula used here is based on the assumption of 
the Bunch-Davies vacuum, the present expansion law is differ- 
ent from the de Sitter one, which may change the expression 
of S H - 
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should caution again that the result of Eq. fl2.8| ) could 
be modified by a careful study of the quantum theory in 
the brane- world inflation and the evolution of cosmologi- 
cal perturbations in the 5-dimensional Einstein equations 

@@. 

W e hav e nu merically solved the b ackg round equations 
( |2.2| ) and (2.3) with the initial value ( |2.6| ), and confirmed 
that analytic estimations are in good agreement with nu- 
merical results (see Fig. 1). 



where and p x are the energy densities of produced par- 
ticles. In the stage where contributions of x are negligible 
in Eqs. (3/2) and ( |3.5| ), making use of the time averaged 
relation, (^(/> 2 )t = (V(4>))t, gives the following approx- 
imate solutions for the scale factor and the background 
inflaton field: 



"' 1 T 



1/3 



<fi(t) = $(t) sinmt, 



(3.6) 



III. BRANE PREHEATING 

Let us investigate a two-field model where the second 
scalar field, \, is coupled to massive inflaton, cj>, with a 
four-leg interaction, \g 2 <t> 2 X 2 '■ 

V(<f>, x ) = \™ 2 <i> 2 + \g 2 4> 2 x 2 + \™\x 2 , (3.1) 

where m x is a mass of the x field. We consider pertur- 
bations of scalar fields on the 3-brane in the flat FRW 
background and implement backreaction effects of pro- 
duced particles at second order JToj ] . 

Decomposing the scalar fields into homogeneous parts 
and fluctuations as 4>(t, x) — > cf)(t)+54>(t, x) and x(i, x) — > 
x(t) + 5x(t> x ); the equation for the homogeneous part of 
inflaton can be written as 



where the amplitude, <&(£), is 
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3/2 



(3.7) 



This evolution is different from the standard results in the 
dust-dominated preheating scenario, a oc t 2 / 3 and $(t) 
t~ 2 1 1 1 . Particle creations in an expanding universe are 
strongly affected by the change of background evolution. 

Introducing a new variable, 8X k — a 3 ^ 2 Sxk, and ne- 
glecting a curvature term in the frequency of SX k which 
is unimportant during reheating, we find that Eq 
is reduced to the Mathieu equation 



d3j) 



£ 

dz 



—SX k + (A k - 2q cos 2z) 6X k = 0, 



(3.8) 



where z = mt and 



. + 3ff0+(TO 2 +.g 2 ( X 2 »0 = O, 



(3.2) 



where (x 2 ) is the expectation value of x 2 , which works 
as backreaction at the final stage of preheating. Neglect- 
ing the contribution of inflaton fluctuations which can be 
negligible in the present model unless rescattering effects 
are taken into account, each Fourier mode of the $Xk 
fluctuation satisfies 



Sxk + 3HSxk 



k 2 



2 X2 



(t) 



S Xk - 0. (3.3) 



The oscillation of the inflaton condensate leads to para- 
metric amplification of the Sxk fluctuation, whose vari- 
ance is defined by 



<<5X 2 }~ / k 2 \Sxl\dk. 



(3.4) 



Note that the equ atio n for x(t) is obtained by substitut- 
ing k = in Eq. (3.3), whose evolution is similar to the 
long wave Sxk modes. 



Assuming that the linear term of p in Eq. (^2) can be 
neglected relative to the quadratic term during the whole 
stage of preheating, the Hubble parameter satisfies 
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(3.5) 
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The efficiency of resonance strongly depends upon the 
initial value of q (— qi). In the standard preheating sce- 
nario where q decreases as q ~ t~ 2 , we require the cou- 
pling, g > 10~ 4 which corresponds to qi > 100, for the 
relevant growth of Sxk fluctuations |l0| . 

In contrast, in the brane-world preheating, q evolves 
more slowly as q ~ t , which indicates that large qi 
much greater than unity is not required for x particle 
productions to occur. This has some analogies with the 
nonminimally coupled x field in the R 2 inflation model 
where slow decrease of the oscillating scalar curvature 
(R oc t^ 1 ) leads to strong preheating EQ]. The ratio 
m/Ms is an important factor to determ ine t he strength of 
resonance. Adopting the value of Eq. (2.S) which comes 
from the CMB constraints, one obtains 



1 



Q ' 



0.36 (10 6 5 ) - 



(3.10) 



The coherent oscillation of inflaton turns on around z ~ 
7r/2, which yields the relation, % ~ (10 6 g) 2 /5. When 
m x = 0, we find numerically that Sxk fluctuations exhibit 
exponential increase for g > 5 x 10~ 6 , corresponding to 

1i > 5 - 
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Let us consider concrete cases with a massless \ field. 
When g — 5 x 10~ 6 , the variance (5x 2 ) increases at the 
initial stage, but its growth stops around mt ~ 35 (see 
Fig. 2). Although \ particles are created while 5\k fluc- 
tuations pass through instability bands with the decrease 
of q due to cosmic expansion, no creation occurs once q 
drops down to 1/4 ~ 1/3, since no instability bands exis t 
there jnj. Actually substituting z = 35 for Eq. ( 3.1C ) 
gives the value, q ~ 0.26, which coincides well with the 
above estimation. 
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FIG. 2. The evolution of the variance, (Sx, 2 ) = {S X 2 )/Mi, 
during preheating for g = 5 X 10 -6 and g = 1 X 10~ 5 
with m x = 0. While resonance soon terminates at the ini- 
tial stage for g — 5 x 1(F 6 , % particles are sufficiently pro- 
duced for g = 1 X 1CP 5 . Inset: {Sx' 2 } vs t for the case of 



2 x Hr d ,m„ 



10m. The x particle whose mass is 



greater than the inflaton mass, m, can be sufficiently created 
for g > 2 x 10" 5 . 

With the increase of g, the duration during which the 
X field stays in instability bands gets longer. In Fig. 2 we 
plot the evolution of (Sx 2 ) for g = 10~ 5 . We find that 
the field reaches a short plateau around 25 < mt < 50 
after some growth in the initial stage. This corresponds 
to the stability region around 1 < q < 2, after which 
(5x 2 ) begins to grow again since the held passes through 
the first instability band around 1/3 < q < 1. In this case 
the backreaction effect of created x particle is marginally 
important at the hnal stage of preheating. 

When g > 10~ 5 , the backreaction effect plays an im- 
portant role to shut off the resonance. While the growth 
rate of Sxk fluctuations gets larger with g being increased, 
the final variance begins to be suppressed, since reso- 
nance becomes less efficient when the g 2 (x 2 ) term in 
Eq. ([T^) grows to the order of to 2 . Numerical simu- 
lations based on the Hartree approximation imply that 
the final variance slowly decreases for g > 10~ 4 . In order 
to estimate the correct size of fluctuations, however, we 
need to include the rescattering effect (i.e., mode-mode 
coupling), which becomes important when particles are 
sufficiently excited pl[ . 

If the x mass is taken into account, this generally works 



as a suppression for an efficient preheating. However, 
as long as the mass effect does not make the x ^e\d 
deviate from instability bands, massive x particles can 
be resonantly amplified. For example x particles with 
mass m x = to are created for g > 2 x 10~ 5 ; and when 
m x — 10m, g > 2 x 10~ 3 (See the inset of Fig. 2). Since 
heavy particles whose masses are a few times greater than 
the inflaton mass are sufficiently produced even when 
g > 10~ 4 , this may play an important role for the suc- 
cess of the GUT scale baryogenesis scenario [ ^T||22| . It 
will be particularly interesting to estimate the baryon 
asymmetry produced in preheating, taking into account 
the decay of x particles as in Ref. ]22] ]. 

The spin 1/2 (and higher) quantum fields can be also 
produced during preheating j2^j24|], although the Pauli's 
exclusion principle restricts the occupation number as 
rife < 1. Consider a massive fermion field tp with a bare 
mass rav,, which is coupled to inflaton with a Yukawa 
interaction, hc^iptp. Then the Dirac equation is written 
in the form: 

(z 7 M V M - TOcff )%l> = 0, (3.11) 
where m e g is the effective mass of fermion, defined by 



m c ff = + hcf>. 



(3.12) 



In an expanding universe with the FRW background, 
fermion production takes place around the region where 
m e g vanishes, which corresponds to the value, <p* — 
—m^p/h p4[ ]. Fermions are periodically created when 
the inflaton passes through 0*, as long as th e re lation 
to,/, < is sat isfied . Substituting Eq. (3.6) with 



Eq. (B/7) for Eq. ( 3.12 ), one finds that the minimum 
effective mass is achieved for mt — 37r/2, yielding the 
upper limit of massive fermion production as 



TO^ 



3/2 



to ~ (10 6 /i)m. 



(3.13) 



This implies that heavy fermions much greater than the 
inflaton mass are easily excited unless h is small as 
h < 10~ 6 . Since the amplitude of inflaton decreases 
slowly relative the standard preheating, this may delay 
the time at which fermion creations are frozen. It will 
be worth investigating how massive fermions produced 
in brane preheating affect the leptogenesis scenario [E5|. 



IV. CONCLUSIONS AND DISCUSSIONS 

We have studied preheating in the 3-brane world in 
the massive chaotic inflationary scenario. While the exis- 
tence of the brane in higher dimensions works to increase 
the amount of inflation via quadratic modifications of en- 
ergy density, its effect in preheating is to make the adi- 
abatic damping of the inflaton slower. Considering the 
X field coupled to inflaton via interaction, \g 2 (j) 2 x 2 i we 
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find that massless x particles exhibit strong amplifica- 
tion even when the coupling g is small as g — 5 x 10 ~ 6 , 
due to the modifications of the background evolution. 
We have also verified that massive scalar and fermionic 
particles heavier than inflaton can be nonperturbatively 
produced, which may play important roles for the baryo- 
and lepto-genesis scenarios. 

Although we have restricted ourselves in the massive 
inflaton case, the self coupling inflationary scenario with 
an effective potential, V{4>) = \^4> A + \g 2 <t> 2 X 2 i has a 
different structure of resonance. I n th is case introducing 
a new variable, SXk = aSxk, Eq. ( |3.3| ) is reduced to the 
following Lame equation in the linear stage of preheating: 



d^ 6xk 



a 2 
A 



i 

71 



sx k = 0, 



(4.1) 



where x = VX^rj and k 2 = k 2 /(\(f> 2 ), with rj = j a _1 dt 
being a conformal time. The cn (x, l/\/2) is the elliptic 
function which is well approximated as cn (a;, 1/V2) ~ 
cos(0.8472x). Since the effect of the adiabatic damping 
of inflaton does not appear in the frequency of SXk, the 
strength of resonance is determined by the ratio, g 2 /X, 
which is the same as the standard scenario discussed in 
Ref. |^6| . Despite this, the slow growth of the scale factor 
makes the redshift of Sxk less efficient due to the relation 
of Sxk — a~ 1 6Xk, which will provide the larger growth 
of the variance, (Sx 2 )- 

The enhancement of field perturbations can lead to the 
growth of metric perturbations |^,^8|, which is partic- 
ularly important for the formation of primordial black 
holes 29 1 . This requires a consistent study of cos- 
mological perturbations including the effect of the 5- 
dimensional bulk, which we leave to future work. 
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